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I. INTRODUCTION 


Modern tactical missiles require increased stand-off ranges 
and need to meet threats from highly maneuverable air targets. 
The high maneuverability of air targets has directed the use of 
defence missiles capable to develop higher lift accelerations and 
more complex control laws. In order to accomplish the 
requirement of large stand-off ranges, propulsion systems using 
air-breathing engines have been studied and developed in recent 
years. The advent of air-breathing engines has naturally led to 
Uae cotsideratvion of BIT missiles in order to minimize the inlet 
angle-of-attack. 

RR nece ss or more complex control laws has introduced the 
IA cation of modern control and estimation theory, since more 
complicated information of the missiles states are needed. BTT 
controlled missiles are generally characterized by increased 
maneuverability and considerable drag reduction over conventional 
Crucemorm, roll» stabasazed STT controlled missiles. Certain 
limitations in technology [Ref. 1] have delayed the development 
of BIT control systems, and consequently any progress in the area 
of BTT autopilots. 

Major technological improvements during the last decade, as 
the availability of advanced digital computers, reopened the 
issue and made BIT control feasible in spite the added complexity 


ОТ Centre] laws@for thecautopidots.: “Tm addition; certain types 


1$ 


of ramjet engines [Ref. 2], which are candidate propulsion 

systems for modern tactical mission requirements of range and 

high altitude [Ref.3], have presented a need for a missile control 
technique to maintain effective inlet flow. This was the main 
reason for given further impetus to the investigation and 
development of BTT control. 

Despite the fact that BIT steering may provide improved 
performance for a missile system, there are still unanswered 
questions concerning stability during homing phase, guidance 
performance, autopilot guidance logic and subsystem requirements. 
All these questions have to be investigated and properly answered 
in order for BIT steering to be considered as a viable control 
method for high performance missiles. 

During the past decade many missile programs [Ref. 3] were 
initiated to improve the capability of steering tactical "egg 
via BTT control with results that have greatly advanced the 
understanding of the various missile subsystems. In the 
autopilot area many different types have been designed and 
developed. All of them force the missile to roll or bank, so that 
the steering maneuver occurs with the missile axis oriented ina 
specific or preferred direction with respect to the incoming 
airstream. This class of autopilots is usually known as POC 
autopilots. 

The main criterion for the seleetion of a particular. SN 
autopilot is based upon the guidance, airframe and propulsion 


system requirements. Generally, missiles with either one or two 
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pwes of symmetry use a POC autopilot which forces the missile 
Ek ООП order to turn as an aircraft. If this motion is 

ER rdinated, then the autopilot is referred to as a CBTT 

ERRO pilot. 

In the guidance area, radome aberation effects for frequency 
guidance are of major concern [Ref. 3] and are being investigated 
in great extent. Also, the interaction between BTT control, 
antenna stabilization and sensor orientation are some of the 
additional concerns that have to be properly addressed. However, 
simplified studies [Ref. 4], which neglect radome effects and 
assume that the missile motion is entirely coordinated, have 
moven that BIT control can provide acceptable performance with 
roll rates that are not excessive for autopilot design. These 
studies were made for a medium range area and long suppression 
me won, and considered both high lift (i.e.; planar) and 
moderate lift (i.e., cruciform) airframe configurations. 

me осе to take full advantage of CBTT control, planar 
airframes have been designed to increase the lifting capability 
in ‘one direction without the weight and drag penalty associated 
with orthogonal lifting surfaces [Ref. 5]. These airframes have 
aerodynamic properties characterized by increased potential to 
ENMGNHNSENUBTT control. 

The present thesis addresses the design and analysis of 
Serete lateral autopilots for application to CBIT missiles. 

The first part reviewed the design procedure of the two 


individual lateral channels as developed by Arrow [Ref. 6]. The 
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design was performed using classical techniques and involved the 
uncoupled yaw and roll channels for the elliptical and cireular 
airframes respectively. The resulted continuous open loop 
designs were analyzed in terms of their transient and frequency 
responses and found to be in accordance with the desired 
requirements specified in [Ref. 6]. Then, applying analog-to- 
digital conversion, the corresponding discrete lateral autop os 
were obtained and analyzed. 

The second part utilized modern control design methods e 
already discussed discrete single-input lateral autopilots. This 
allowed comparison with the preceding classical design, and more 
importantly established a technique to extend some of the results 
to the more general multivariable case. At first, assuming 
availability of all states for feedback purposes, application of 
the Ackermann formula led to a discrete state-feedback designed 
autopilot. Next, since the state vector of the state-feedback 
model is not usually accessible to direct measurement, an estimator 
was introduced as an additional dynamic design in order to 
implement control to the original system. The state-feedback and 
estimator designed autopilots were analyzed for both lateral 
channels and found to have satisfactory responses. 

Finally, coupling the discrete pitch and roll ekannali 
autopilots, the state-feedback and estimator designs were 
obtained and proved to be robust. 

The analysis in all the above cases was performed using the 


existed at Naval Postgraduate School OPTSYS and ORACLS Forte 
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CO On Encinas and discrete systems respectively. 
Additionally, the last part that dealt with the coupled pitch and 
men channel autopilot utilized the POPLAR design program 


developed by Gordon [Ref. 7]. 
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II. CLASSICAL DESIGN AND ANALYSIS OF LINEAR UNCOUPLED 
LATERAL AUTO ripen 
A. GENERAL 
The initial phase in the design of lateral CBTT autopilots 
involved the design and analysis of the individual uncoupled 
lateral channels (i.e., yaw and roll) with prescribed 
relationships between speeds-of-response. These relationships 
when coupled with the corresponding ones of the longitudinal 
uncoupled channel (i.e., pitch) would meet the requirements of 
the overall-CBI autopilot. 
The uncoupled autopilot design method was classical and used 
a combination of frequency response and root locus techniques. 
Utilization of this particular design method Wed tomer 
achievement of practical bandwidths (i.e., sufficient high 
frequency attenuation), and.in turn provided the range of required 
missile body angular rates and control motions. In addition, the 
resulting design minimized the influence of aerodynamic 
variations on desired responses. The application of the 
uncoupled channels to the whole CBTT autopilot was accomplished 
by an appropriate choice of the relative time constants of the 
individual channels. In order to achieve the desired maneuver 
plane acceleration the roll channel was designed to have a time 
constant of 0.5 seconds. The yaw uncoupled channel, which 
follows the roll motion to produce the required coordination 


(i.e., minimization of sideslip angle), is designed to have a 
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more rapid response with a time constant of 0.39 seconds for the 
circular airframe. The detailed requirements for the classical 
design of the uncoupled yaw and roll channel autopilots are 
presented in Appendix A. 

A fixed flight condition (i.e., constant Mach number and 
altitude) was selected for this preliminary performance study. 
Fixed flight conditions are typically used in autopilot designs 
to identify and cure critical areas of concern. When autopilot 
requirements are satisfied at fixed flight conditions, then areas 
of concern caused by varying them are addressed. The selected 
flight condition of 60000 feet altitude and Mach number 3.95 
provided sufficient dynamic pressure, so that the missile 
maneuvers resulted in large enough angles-of-attack to exercise 
Sis Up control) Aerodynamic data for this particular flight 
condition are provided in Appendix A. 

The aerodynamic models developed for stability studies in the 


frequency domain were linearized about a trim angle-of-attack for 


both lateral channels. The following three assumptions were 
made: 
1. The plane ^g - Ze of Figure 2.1 was the maneuver plane. 


2. The missile was trimmed in pitch (i.e., M yz O, at fixed 
values of «X , q, апа ©. 
Rather than use the assumption that the missile roll rate is 
P or mately Zero as rt 1s"normally done for the roll stabilized 


STT control, the following assumption was made for BTT: 
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se missile roll rate was constant. 
Linearized aerodynamic derivatives are given in Appendix C. 

In this chapter the analysis of both continuous and discrete 
uncoupled lateral channels of a BIT autopilot was based on the 
transient and frequency responses of maneuver plane accelerations, 
body angular rates and tail incidence angles. A general block 
Ian of a BIT autopilot with all its channels is shown in 
Figure 2.2. Inertial acceleration command were applied in polar 
coordinates (i.e., magnitude of the command uM applied to pitch 
and the direction Œ. erol autopilot). The yaw autopilot was 
IEEE G the roll autopilot in order to minimize the yaw and 
roll motions. Achieved maneuver plane accelerations in 
rectangular coordinates (i.e., Ng and Ny ) were determined by 
resolving achieved body-fixed accelerations (i.e., Ng and il, ) 
through missile roll rate (i.e., Euler angles 9 and W were 


assumed to be sufficiently small). 


B. AIRFRAME CONFIGURATIONS 

The two airframe configurations studied in this work were 
taken from [Ref. 6] and are shown in Figure 2.3 and Figure 2.4. 
Although the configuration in Figure 2.3 reveals a body of 
Circular cross section and that of Figure 2.4 an elliptical one, 
both airframes have the same cross sectional area distribution. 
liu corre, theweircular cross Sectional body has a closure 
ratio Abase/Amax of 0.69 with Amax occurring at 68% missile body, 


whereas the elliptical airframe has a 3:1 cross section. 
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EKI conrrguratrons are tail-controlled using four 
ERR mica] control Surfaces which are located flush with the body 
base with a +30° dihedral. Їп the case of the elliptical body, 
the hinge line was skewed such that a 10° control deflection 
measured at the body-tail juncture had a resultant 7.04° surface 
deflection. Thus, the aerodynamic control effectiveness in terms 
of deflection measured at the body-tail is lower for the 
elliptical airframe although it is nearly the same in terms of 
resultant surface deflection. 

The total span of the mono-wings is the same for each 
Generation; which results in larger wing area for the circular 
airframe. The wing area and span for the circular airframe were 
emeeen as typical of current maneuvering missiles. The wing for 
pe ptical concept was determined Dy projecting the 
elliptical body on the circular body-wing planform. The resultant 
when exposed wing planform became the wing for the elliptical 
body. 

Comparison of the elliptical airframe with the corresponding 
circular indicates the following: 

1. About 30% more normal force that is nearly independent of 
angle-of-attack can be achieved at supersonic speeds. 

FE ес od VPongitudinal stability parameter COMME more 
positive, and with more pronounced nonlinearities in 
pitching moment at subsonic speeds. 

3. Levels of directional stability are increased and more 


compatible with levels of longitudinal stability. 


sr 


4. More yaw control is available although suitable locations 
for tails on the body are more limited because of the 
geometry of the elliptical airframe. 

The two airframe configurations were sized to provide 
realistic geometric and mass properties. The details are 


presented in Appendix D. 


C. UNCOUPLED YAW CHANNEL AUTOPILOT FOR ELLIPTICAL AIRFRAME 

The purpose of the uncoupled yaw channel autopilot of a CBTT 
missile is to minimize the sideslip angle (8), or provide 
coordinate motion between the yaw and roll channels. The easiest 
way to accomplish this is by designing the uncoupled yaw channel 
(i.e., roll and pitch dynamic effects neglected) as a regulator 
(i.e., no guidance command and with rate and acceleration 
feedback) to help minimize the sideslip angle. 

A block diagram of the uncoupled yaw channel is shown in 
Figure 2.5. In this diagram both the aerodynamic model and yaw 
control law are involved. The normal acceleration (Hy is not 
used to command the CBTT autopilot. Instead, it is used TY 
design and analysis of the uncoupled channel. The command used by 
the coupled system is shown in dashed lines and is a yaw angular 
rate command (ro). The yaw control law shown in Figure 2.6 [Ref. 
6] is governed by missile body angular rate (r) and yaw nonm Mi 
acceleration Lei At the flight condition of interest (i.e., 
60 kft altitude and Mach number 3.95) tne yan contro law 


determines the required command (с ) to an actuator which is 
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approximated as a first-order lag at 30 Hz. The rate compensator 
computes the high frequency attenuation and is used to minimize 
aerodynamic variations on the quality of the regulator. On the 
other hand, the acceleration compensator measures the acceleration 
bandwidth via the time constant of the acceleration response of My - 
er ansterkunetions of. Aérodynamic Model | 


The aerodynamic transfer functions of the uncoupled yaw 


channel autopilot are: 
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(573) 454 


177 


oo (II.C.1-6) 


> 
K=: ——_— (II C ME 
ү 
Substituting the values of aerodynamic data (Table VIT), 
linearized aerodynamic derivatives (Table VIII), and geometric 


and mass properties (Table IX), equations (II.C.1-3) through 


(II.C.1-7) become: 


A = -0.1351 (II SC. eo 
B = 0.0436 Quo |) 
GE 1 2g (II CO 
E = 34.5495 (Lt. ZER 
[SAS (II.C.1-12) 


Introducing the above equations (II.C.1-8) through 
(II.C.1-12) to (II.C.1-1) and (II.C.1-2), the aerod yrange 
transfer functions of the uncoupled yaw channel for the circular 
airframe and zero angle-of-attack can be obtained. 

a. Transfer Function of Yaw Angular Rate: 


= 
= m S t1 
= 005 (Gai) (II.C. SES 


(a) (as) 





C4 


DE Transfer Function of Yaw Normal Acceleration: 


e (¡As (sn) 
hy ade Gere ко s 


E 1 l | 
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Heameamemne the above transfer functions in terms of the 


diable pairs Ca ea) and (ny , oy) respectively, they become: 


S$ 40.0494 25 -19,1S15 t= - 36,66 ys -4, 74225, (II.C. 1-15) 


2 
Ny S 0.0494 п ү5-18281$п ү: 2.4625 5, - 0.0018 oy $- 524, 6801 oy imc 1-16) 


Applying inverse Laplace transformation, the following 


set of linear differential equations oan be obtained: 


MA -15.2915 7s - 36.66 by - 3.7422 50 (II.C.1-17) 


Ny 0.0494 ny -18, 2515 ny = 2.4625 by - 0.0018 dy - 524. 6Єбс18, (11.c.1-18) 


Both equations (II.C.1-17) and (II.C.1-18) form a second- 
order system of linear differential equations in which the 
forcing function involves derivative terms. Using rules of 


state-space representation [Ref. 8] the following equations are 


obtained: 

Kz X,- 56.66 Sy mt Lo 
xy 3 18,2815 X,-0.0494 X,-=0, 9312 cy (II.C.1-20) 
ee se ew. (II.C.1-21) 
Z * 18.4815 Z -0.0494 2, - 280.6558 e. (II.C.1-22) 
^ Quien 3) 
€ z X,- 56.66 су Gneo) 


Ras 5, + 2.4625 7 (HE. 1-25) 


il 


E, 


2. Equations of Yaw Control Lav onee eon 
a. Acceleration Compensator Equation 

The acceleration compensator equation obtained from 

E es. 6 is: 
= 031746 Ghisa ) (II.C. 25 
OAS ti Le 

Rearranging and applying inverse Laplace transforma- 
tion, (II.C.2-1) turns into the following linear different mem 
equat Tom: 
y==5y PST EE We (II 20 

Substituting equation (II.C.1-25) into (IE.C. 2 
the last becomes: 
y=1,597732 -5y+ 3,7554-cy - 1.5713 My, (II. Cc. 2868 

b. Rate Compensator Equation 

The rate compensator equation also obtained from 
Figure 2.6 D | 
" 4 35 4.95 Gest) (II.C.2-4) 
M. Å Cy TE) 

Rearranging and applying inverse Laplace transforma- 
tion to (II. CRER st turns into: 


0.485 у+0.4@5 2 +4855у+455 (ттс ШЕ 


dd 
Substituting equations (II.C.1-23), (II.C.1-20]m P 
(II.C.2-3) into the above, it becomes: 


ey = $85 X, FOES XK tC 4TZ 4 24457215, $124 oy - 0.7747 ny (ICAA 
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C. Actuator Equation 


The actuator equation obtained from Figure 2.5 is: 





— ] x 
E - E = M gu С. 2—7} 
— stl 
188.4 


Rearranging and applying inverse Laplace transforma- 


kwen, the last equation turns into: 


3. 188 4è, -18848, fer .c.2-8) 
É 


3. Design Approach and Analysis of Continuous System 


Utilizing state-space representation, the equations 
EN 1-19) through (II.C.1-22), (EII.C.2-2) (II.C.2-6) and 
(II.C.2-8) ean be modeled in a seventh-order system of the form 
Meee. The continuous plant system and input matrixes F and G 


are shown in Table I, and the state vector is: 


di (TMC. 3-1) 


5; 


doere the state variables are: 
X1, X5: yaw angular rates 
Z1, Z9: yaw normal accelerations 
y : output of acceleration compensator 
Sy : input command in the actuator 


yaw tail incidence 
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MANYOK Lhe pole-zero and time and frequency response plots are 
ОЕ пед. 

The pole-zero plot of Figure 2.7 indicates that the 


continuous open loop system is stable, since the s-plane poles 


ake. 

SL STE (II.C.3-2) 
Sz = - 6.1042 tj iC. C $9€ (II.C.3-3) 
53 =-6.1042 - j 10.6376 (II.C.3-H) 
$4 7 -O.C2174T$ * j 08914 Eno 3-5) 

ET T - i| OG 

SE -0.0217475-J 4.0891 ke) 
DEL 0776 12, 729 29 Lon mt 
57 :-4.9296 - 12.949928 CIS =0 ) 


The time response plots of the yaw normal acceleration, angular 
rate and tail incidence are shown in Figures 2.8 through 2.10. 
In particular the yaw normal acceleration time response plot has 
a 0.39 seconds time constant, 7% overshoot and a steady-state 
error of 0.018. These results are in accordance with the 
requirements referred in Appendix A, that is a time constant of 
O.4 seconds, overshoot less than 10% and steady-state error not 
necessarily equal to zero. All the above three time response 
plots are identical with those presented in [Ref. 6]. 

Figures 2.11 through 2.16 show the frequency response 


plots of the yaw normal acceleration, angular rate and tail 
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incidence, from which the phase crossover frequencies and gain 
margins of Table II can be obtained. The positive gain and phase 
margins of the open loop system ensure the relative stability of 
tne closed loop (controlled) system. 
4. Design Approach and Analysis of Discrete System 

Utilizing analog-to-digital conversion by the aid of 
ORACLS program and for a sample period of 0.0125 seconds, a 
seventh-order discrete system of the form x(k+1)=Ax(k)+Bu(k) is 
obtained. The discrete pant system and input matrices A ann 
are shown in Table III. 

The pole-zero plot of Figure 2.17 indicates that the 
discrete open loop system is also stable, since the z-plane 


poles are: 


Z,= 0.113076 (II.C.4-1) 
Zz= C.91E354+|0.1228€2 (II. ON EE 
Zz = CIE SSA OE | (II. Cleo 
Z4: 0996422 4 JO.OS10768 (II.C.4-4) 
Zs = O.996 422 -] COSIOTSS (II.C.4-5) 
. . ONA EE 
D = OS JONE (II GE 
< L ese ©) Se = 
Z7 *0. 96 556S - 0.036 1346 (II.C. 4-7) 


The time response plots of the yaw normal acceleration, 
angular rate and tail incidence for the discrete uncoupled nu 
channel are presented in Figures 2.18 through 2.20. A close 
observation of these plots indicates that they are identical with 
those of the continuous classical system found in the previous 


section. 
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Figure 2.10 Yaw Tail Incidence vs Time, Uncoupled Yaw Channel 
Autopilot; Classical Design; Continuous Open 
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Figure 2.11 Yaw Normal Acceleration-Gain vs Frequency; Uncoupled 
Yaw Channel Autopilot; Classical Design; Continuous 
Open Loop System; Elliptical Airframe 
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Figure 2.15 Yaw Tail Incidence-Gain vs Frequency: Uncoupled 
Yaw Channel Autopilot, Classical Design: 
Continuous Open Loop System: Elliptical 
Airframe 
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D. UNCOUPLED ROLL CHANNEL AUTOPELCT FEE EEROR eo 

The uncoupled roll channel autopilot of a CBITT missi lem 
commanded to roll the missile so as to put the preferred maneuver 
direction in the direction of the guidance acceleration command. 
The desired maneuver plane acceleration should be attained as 
rapidly as the achieved body-fixed pitch acceleration. To 
accomplish this, the uncoupled roll channel autopilot (i.e., yaw 
and roll dynamic effects neglected) was designed to have the roll 
angle time constant equal to the time constant of the normal 
acceleration achieved by the uncoupled pitch channel autopilot. 

A block diagram of the uncoupled roll channel is shown in 
Figure 2.21. In this diagram both the aerodynamic model 3mm 
control law are involved. The roll control law shown in EM 
2.22 is commanded by roll angle (0,) and governed by roll angular 
rate (p) and roll angle (0). 

The design and analysis of the uncoupled roll cnar i 
autopilot was performed in this section for the stable at zero 
angle-of-attack circular airframe. 

1. Transfer Functions of Aerodynamic Model 

The aerodynamic transfer functions of the uncoupled roll 


channel autopilot obtained from Figure 2.21 are: 


a. Transfer Function of Roll Angular Rate: 
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Figure 2.17 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 
Classical Design; Discrete Open Loop system; 
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Figure 2.19 Yaw Angular Rate vs Time; Uncoupled Yaw 
Channel Autopilot; Classical Design: 
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b. Transfer Function E 
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а (II.D MEN 
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Substituting the values of aerodynamic data (Table VII), 
linearized aerodynamic derivatives (Table VIII) and geometric and 
mass properties (Table IX), equation (II.D.1-1) becomes 

ps = 8.03462 Cp (II.D. PS 

Applying inverse Laplace transformation to equations 

(II.D.1-2) and (II.D.1-3), the following set of linear 


differential equations is obtained: 


p = 8.03462 Sg (Il Dee 


DE Dp (Il De 5) 
2. Equations of Roll Control Law and Actuator 
a. Roll Angle Compensator Equation 
The roll angle compensator equation obtained from 


Figure 2.22 99m 


C 
AE «<< Sea (11 EE 


с 
Rearranging and applying inverse Laplace 
transformation (IT.D.2-1) turns into the following lined 
differential equation: 
X = -17.69- 8X + 17.69, ae c 
b. Rate Compensator Equation 


The rate compensator equation also obtained from 
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aan nd applying inverse Laplace transforma- 
КООП to (IT.D.2-3), it turns into: 
Y + 5Y = 0.0502X - 0.0502p + 0.75325X - 0.75325p gan 224) 
EE a ecquations (IT.D.1-4) and (11.D.2-2) into 
(II.D.2-4), the last becomes: 
Y2-0.7532p-0. 88352440.35165X-5Y-0. 403340, +0.883524, DOS) 
cC.  Pseudo-Differential Equation 


ООШ О МОЛО а Геге ла eadauation obtained from Figure 


ame? is: 
0. C S 
En "o7 bE 2-6 
боозу 
a 
Rearranging and applying inverse Laplace 
ea formation to (II.D.2-6), it turns into: 
X, = -6X, + 0.078198p (II.D.2-7) 


ОШ ng equation (TI.D.1-4) into (II.D.2-7), the 
last becomes: 
X. = - 6X, + 0.628291 cp Gu iS) 
ata tion OF AHetuator Compensator 
The equation of actuator compensator obtained from 


ШКОЛ ге 2.22 is: 


> 
no + 
5 
i c 


1 


C Y-X,) Gun D 2-9) 


Vi 


Rearranging and applying inverse Laplace transforma- 
Baao (II.D.289),..dt turns into: 


23 =-15 Cp, +0.13636Y-0.13636X,+15Y-15X, (II.D.2-10) 
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Substituting equations (II.D.2-5) and (11.D.2-8) ЧИЕ 


(II.D.2-10), the last becomes: 
ER --CJAC27I p -0.RC4S q eC. CAT 3 X4: 3182 Y -IEISIE X “15 eg, - Q [ICE fc tE. 
E 
(III. DENE 


e. Actuator Equation 


The actuator equation obtained from Figure 2.2 ER 





= SAZ 
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[EE 4 


Rearranging and applying inverse Laplace 


transformation, еш ARNT NN: 


E 10795. 326, -188.1 Op mme 
C 


C4- 
T 


3. Design Approach and Analysis of Continuous System 


Utilizing state-space representation, the equations 
(II.D.1-4), (II.D.1-5), (II.D.2-2), (II.D.2: (IT. DOM 
(II.D.2-11), and (II.D.2-13) can be modeled in a seventh-order 
System of the form x=Fx+Gu. The continuous plant system and 


input matrices F and G are shown in Table IV, and the stan 


vector is: 
f 


ACK) = 


ER SOUS 


(II.D = 


QA ^ 
A 


Wa 
A 
TE 
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where the state variables are: 

roll angular rate 

roll angle 

output of roll angle compensator 


output or rate compensator 


PS Hq bi «© ‘0 


output of pseudo-differentiator 


input command in the actuator 


C1 
A 
(Ga 


roll tail incidence 


O4 
= 


Ementa che OPTSYS program, using an input step 
Amon which represents "| zee command" at zero trim angle-of- 
attack, the pole-zero and time and frequency response plots are 
obtained. 

Mie. poleszero plot oi Figure 2.23 indicates that the 


continuous open loop system is stable, since the s-plane poles 


are: 
3, m. 785 (roll angular rate) QU Do SA) 
>> = -9.25097+j28.4098 (roll angle) (II.D. 3-3) 
Sz = -9.25097-j28.4098 (output of roll angle compensator) 
CHIN 
S4 = «2.46608+j2.71152 (output of rate compensator network) 
(IT.D. 3-5) 
AE 2. 16000-92.71152 (output of pseudo-differentiator) 
(II.D.3-6) 
> = -8.98209 (input command in the actuator) mer) 
57 2:165 19032 (rol lea) i memdence ) (I.D. 3-85) 
The time response plots of the roll angle, angular rate 
Emu iIncagdence are shown in Figure 2.24 through 2.26. In 
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Ia FU HOI angle time response plot has a 0.55 seconds 
time constant, 3% overshoot and a steady-state equal to zero. 
These results are in accordance with the requirements referred in 
Appendix A, that is a time constant of 0.5 seconds, overshoot 
less than 10% and zero steady-state roll angle error. All the 
above three time response plots are identical with those 
presented in [Ref. 6]. 

Figures 2.27 through 2.32 show the frequency response 
plots of the roll angle, angular rate and tail incidence, from 
which the phase crossover frequencies and gain margins of Table 
V can be obtained. The positive gain and phase margins of the 
open loop system ensure the relative stability of the closed loop 


(controlled) system. 
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PHASE CROSSOVER GAIN 
FREQUENCY (rad/sec) MARGIN (db) 
ROLL ANGLE () 10.0461 60.2165 
ROLL ANGULAR RATE (p) 10.0461 60.2165 


ROLL TAIL INCIDENCE (de ) De A 
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4. Design Approach and Analysis of Discrete System 


Utilizing analog-to-digital conversion by the aid of 
RES program and for a sample period of 0.0125 seconds, a 
seventh-order discrete system of the form x(K+1)=Ax(k)+ Bulk) is 
obtained. The discrete plant system and input matrices A and B 
are shown in Table VI. 

The pole-zero plot of Figure 2.33 indicates that the 


discrete open loop system is also stable, since the z-plane poles 


are: 
=, = 0.0992805 (roll angular rate) GET E e 
Z, = 0.835216+j0.309736 (roll angle) (II.D.4-2) 
Zz = 0.835216-j0.309736 (output of roll angle compensator) 
(II.D. 4-3) 
Za = 0.969087+J0.0328588 (output of rate compensator network) 
E TD AI) 
Zs = 0.969081-j0.0328588 (output of pseudo-differentiator) 
(II.D.4-5) 
Ze = 0.893797 (input command in the actuator) (II.D.4-6) 
Z, = 0.937075 (roll tail Mce) (II.D.4-7) 


Ine time response plets of the roll angle, angular rate 
and Gail incidence for the diserete uncoupled roll channel are 
presented in Figures 2.34 through 2.36. A close observation of 
Bae plok. indicates that theyeare identical with those of the 


Continuous classical system found in the previous section. 
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Autopilot; Classical Design; Discrete Open Loop 
System; Circular Airframe 
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DI MODERN CONTROL DESIGN AND ANALYSIS OF LINEAR 
UNCOUPLED LATERAL AUTOPILOTS 
A. GENERAL 

The task of this chapter is the design and analysis of the 
Same discrete uncoupled lateral channel autopilots discussed in 
the previous chapter, using different techniques which are based 
On modern control formulation. The difference in the two 
BHEDDOSches is entirely in the design method since the end result, 
CET of difference equations providing control, is identical. 

Modern control theory is contrasted with the classical 
control theory in that the former is applicable to multi-input- 
multi-output systems, which may be linear or nonlinear time- 
invariant or time-varying, while the latter is applicable only to 
linear time-invariant single-input-single-output systems. Also, 
modern control theory is essentially a time-domain approach, 
While the conventional classical control theory is a complex 
frequency-domain approach. 

System design in classical control theory is based on trial- 
Werner procedures which, in general, will not yield optimal 
se nos stemmer system design in modern control theory, on the 
other hand, enables the design of optimal control systems of 
great complexity and good accuracy with respect to given 
КООШО? Се indexes. In addition, design in modern control 


ШОГУ Can be carried out for a class of inputs instead of 
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specific input function, such as the impulse, step or sinuso 
functions and can also include initial conditions. 

One of the most attractive features of modern control design 
method is that the procedure consists of two independent steps. 
One step assumes that all the system states are available for 
feedback purposes. In general, even if this is not a practical 
enough assumption since it needs a large number of sensors, it is 
usually adopted in order to accomplish the first design step; 
namely the control-law. The remaining step is the design of an 
estimator which estimates the entire state vector, given 
measurements of portion of the state provided by the system output 
equation. The final control algorithm consists of the control- 
law and estimator combined, where the control-law calculations 
are based on the estimated states rather than the actual states. 
This substitution is reasonable and the combined design can give 
closed loop characteristics which are unchanged from those 


assumed in designing the control-law and estimator separately. 


B. DISCRETENSTANES-EEEDBACK DES EN 
Considering the following discrete control system: 

X(K+1) = Ax(k) + Bulk) (ттт: | 

MO labo (ke) (III.B.1-2) 
the control-law design is also referred to as state-feedback 
design since it is simply the feedback of a linear combination: 
all the system states, that is: 

u(k) = Fx(k) (ТТТ.В.1-3) 


where F: control-law gain vector 
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Thus, the characteristic equation of the controlled (closed loop) 
system is: 
det (ZI-A+BF) z 0 CAT dee ) 
The discrete state-feedback design, providing that the system is 
controllable, consists then of finding the control-law gain 
Wélekor E so that the roots of (III.B.1-4) are in desirable 
ieeations. 
iWibresram logie for computing the control-law gain vector via 
the Ackermann's formula is given in Appendix E [Ref. 9]. 
paang this control algorithm a Fortran program was written 
(Appendix F) which has as inputs the sample period, the discrete 
plant system and compredl input matrices, the s-planeepoles and 
provides as output the control-law gain vector. 
КОО еен реа Yaw Channel for Elliptical Airframe 
a. Control-Law Gain Vector 
Executing the Ackermann Fortran program described in 
Appendix” P with inputs: 
(1) Sample period of 0.0125 seconds 


(2) Discrete plant system and control input matrices A and B 
of Table III o a 


(3) S-plane poles defined in equations (II.C.3-2) through 
(ТФ®С.3-8) 


the following control-law gain vector for the elliptical airframe 
of the uncoupled yaw channel is obtained: 


F=[-1.0195 -0.109 -0.3492 0,0393 -0.6152 32.4775 -31.0868] 
IIS ES) 
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b. Design Approach and Analysis 
The discrete state-feedback designed yaw autopilot 
can be found by introducing the control-law gain vector of 
(II.B.1-5) into the original system, 
The pole-zero plot of Figure 3.1 indicates that the 
discrete closed lop system is stable, since the z-plane closed 


loop poles are: 


N 
T 


0.46461 (yaw angular rate) (Illi . Beier 
“, = 0.898888+j0. 137039 (III BO 


Zz = 0.898388-j0.137039 (yaw normal acceleration) QRR E WE 8 


Z4 = 0.998423+j0.051078 (III.B.1-9) 
5 = 0.998423-j0.0510789 (output of acceleration compensator 
network) (ILL. BEN 
Z, = 0.962974+j0.0358712 (input command in the actuator) 
(III.B.1-11) 
Zy = 0.962974-j0.0358712 (yaw tail incidence) (III BRENN 


The time response plots of the yaw normal 
acceleration, angular rate and tail incidence are presented in 
Figures 3.2 through 3.4. A close observation of the above pole- 
zero and time response plots for the discrete yaw state-feedback 
design indicates that they are identical with those of the 
discrete classical design found in the previous chapter. 

c.  Simplified Design 

The discrete state-feedback designed autopilot of the 
previous section can be simplified by reducing the returning gain 
loops. This can be accomplished by placing Zeros into apr... 


elements of the control-law gain vector: 
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F=[-1.0195 -0.109 -0.3492 O -0.6152 32.4775 -31.0868] 
9 Klees 13) 


The poleszero and time response plots of the 
resulting simplified state-feedback yaw autopilot, shown in 
Figures 3.5 through 3.8, do not present significant differences 
from the corresponding plots of the discrete classical design 
apart from the overshoot which was slightly increased. The 
discrete closed loop system is again stable, since the z-plane 


closed loop poles are: 


Z, = 0.467984 A ISI) 
Z, = 0.891328+j0.126019 | (Тт E15 ) 
Ze = 0.891328-j0.126019 Cine - 16 ) 
Z4, = 0.998422+ j0.0510788 (ТТ ., 1—17) 
Zs = 0.998422-j0.0510788 (Te. 18) 
Z- = 0.969075+j0.0429695 (ТИР. 619) 
Zr = 0.969075-j0.0429695 (TI Tyee 1 =20 ) 
se incouplea Roll Channel for Circular Airframe 
a. Control-Law Gain Vector 


Following the same procedure as in the yaw channel 
case apart from the use of discrete matrices A and B from Table 
III and s-plane poles from (II.D.3-2) through (II.D.3-8), the 
control-law gain vector for the circular airframe roll channel 
was found to be: 


F=[0 0.0004 0 0.0178 -0.0179 -1.499 0.0023] (TEB 2>1) 
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Figure 3.1 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 
State-Feedback Design; Discrete Closed Loop System; 
Elliptical Airframe 
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Figure 3.2 Yaw Normal Acceleration vs Time, Uncoupled Yaw 
Channel Autopilot, State-Feedback Design, Discrete 
Closed Loop System: Elliptical Airframe 
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Figure 3.3 Yaw Angular Rate vs Time, Uncoupled Yaw Channel 
Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Elliptical Airframe 
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Figure 3.4 Yaw Tail Incidence vs Time: Uncoupled Yaw Channel 
Autopilot, State-Feedback Design: Discrete 
Closed Loop System: Elliptical Airframe 
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Figure 3.5 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 
Simplified State-Feedback Design; Discrete Closed 
Loop System; Elliptical EEUU mE 
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Figure 3.6 Yaw Normal Acceleration vs Time; 
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Figure 3.7 Yaw Angular Rate vs Time; Uncoupled Yaw 
Channel Autopilot; Simplified State- 
Feedback Design; Discrete Closed Loop 
System; Elliptical Airframe 
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b. Design Approach and Analysis 

The discrete state-feedback designed roll autopilot 
can be found by introducing the control-layn gain vector van 
(III.B.2-1) into the original system. 

The pole-zero plot of Figure 3.9 indicates that the 
discrete closed loop system is stable, since the z-plane closed 
loop poles are: 

Zi 80409899714 (noli angular rate) (IITB 92220 
= 0.832909+j0.315088 (roll angle) (Ттт. е 


Zz = 0.832909-j0.315088 (output of roll angle compensator) 
(III.B.2-4) 


Za = 0.969236+j0.0333773 (output of rate compensator network) 
(III.B.2-5) 


la = 0.969236-j0.0333773 (output of pseudo-differentiator) 
(III 220) 


Ze = 0.896583 (input command in the actuator) L LTE S) 
Zr = 0.938233 (roll tail incidence) (IIL. Be eee 

The time response plots of the roll angle, angular 
rate and tail incidence are presented in Figures 3.10 Unpoman 
3.12. A close observation of the above pole-zero and time 
response plots for the discrete roll state-feedback design 
indicates that they are identical with those of the discrete 
classical design found in the previous chapter. 

c. Simplified Design 

The state-feedback roll designed autopilot can be 

simplified by reducing the returning gain loops as follow: 


F=[0 0 0 0.178 -0.0179 -1.499 0] (ттт. ЕЕ 
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Ba E-7ero and time response plots of the 
resulting simplified state-feedback roll autopilot, shown in 
Figures 3.13 through 3.16, do not present significant differences 
rom the corresponding plots of the discrete classical design. 
The discrete closed loop system is again stable, since the z- 


plane closed loop poles are: 


RO 0992715 (FIT ETO) 
Z, = 0.836632+j0.305719 (III.B.2-11) 
EE 0059566322 j0.305719 (ТЕТЕ 12) 
Z4 = 0.967647+j0.0334948 (Lil see2- 13) 
Z5 = 0.967647-j0.0334948 (III.B.2-14) 
ZEE, 695937 (TE, 2-15) 
Z Ez 10.938256 (TA 10) 


eee ES CRETE ESTIMATOR DESIGN 

The state-feedback design discussed in the last section 
assumed that all system eee were available for feedback 
purposes. Since the state vector is not always accessible to 
direct measurement, an estimator is going to De Introducer in 
this section as an additional dynamic design in order to 
implement control to the original system. The estimator design 
method consists mainly of determining algorithms which will 
reconstruct all the states, given measurements of a portion of 


them. 
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Figure 3.9 Pole-Zero Plot; Uncoupled Roll Channel Autopilot: 
State-Feedback Design; Discrete Closed Loop System; 
Circular Airframe 
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Figure 3.10 Roll Angle vs Time, Uncoupled Roll Channel Autopilot, 
Shane NcedbackeDoe neng eDiserete Closed Loop System; 
Circular Airframe 
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Figure 3.12 Roll Tail Incidence vs Time; Uncoupled Roll Channel 
Autopilot; State-Feedback Design; Discrete Closed 
Loop Syswan; Cireubar Airframe 
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Figure 3.13 Pole-Zero Plot; Uncoupled Roll Channel teg 
Simplified State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Considering the same as in the state-feedback case discrete 
control system, a prediction estimator defined by the following 
equation 18*Iintreghecar 

KAKAI) = AK(K) + Bulk) + KØGE) - HX(k)] (III.C.1-1) 
where x: estimate state vector 

K: estimate gain vector 

In this closed loop estimator, shown in Figure Ci ci 
difference between the measured and estimated output is fed back 
and the model is constantly corrected with this error signal 
which is defined as X=x-X. The difference equation describing 
the behavior of the error is obtained by subtracting equation 
from the actual plant output equation (11.8. 1=2): 

x(K+1) = [A - KH] x(k) (Ттт. ШЕ 
Thus the characteristic equation of the controlled (closed loop) 
system is: 

det (zi. = А a 0 (TIL ICO 


The discrete estimator design, providing that the system is 


observable, consists then of finding the estimator gain vector 


{ч 


so that the roots of (III.C.1-3) are at debirable locataom 
The estimator gain vector can be obtained again, as in the 
case of the state-feedback design, by application of the 
Ackermann program of Appendix F with inputs the sample period, 
the transposes of the discrete plant and output matrices, and 


faster s-plane poles from those of the continuous system. 
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1. Uncoupled Yaw Channel fore lp Dc NE 
a. Estimator Gain Vector 
Executing the Ackermann Fortran program (Appendix F) 
wich inputs: 
(1) Sample period of 0.0125 seconds 
(2) The transposes of the discrete yaw system plant and output 


matrices, that is AÍ and H,!, where: 


(III, CREED 


T 
bi 
Oonddao — 


(3) S-plane poles slightly faster than those of the COM e 


open loop system, that is with more negative real parts: 


ag Me (ттт ee 
S, = -6.1142+j10.6396 . (III.C NN 
Se = -6.1142-510.6396 (ILE. Cie 
54, = 20.022 7 0,809909 (III.C.1-8) 
Ss = -0.0227475-j4.08919 (III.C.1-9) 
5& = -2.9396+j2.99929 (ттт. ЕЕ 
57 = -2.9396-j2.99929 (III Gelb 


the transpose of the estimator gain vector for the elle 
airframe of the uncoupled yaw channel is calculated as output, 


from whichaegkhe following K is obtained: 
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-0.0034 
0.0008 
-0.0003 
K = -0.0177 RIES ) 
080005 
-0.0017 
-0.0240 
b. Design Approach and Analysis 
The discrete estimator designed yaw autopilot can be 
found by introducing the estimator gain vector of (III.C.1-12) 
mu the original system. 
The pole-zero plot of Figure 3.18 indicates that the 
discrete closed loop system is stable, since the z-plane closed 


loop poles were found to be: 


E 


| 0.117115 (yaw angular rate) (a NO . 1-13) 


Z> = 0.918174+j0.122926 Gre D 


Zz = 0.918174-j0.122926 (yaw normal acceleration) (Tas GS) 
<a 


= 0.998399+j0.0510958 (TITI 601-10) 
Zs = 0.998399-j0.0510958 (output of acceleration compensator 
network) (Тт 1—17) 
25 


0.963249+j0.0361306 (input command in the actuator) 
GOL CCELI 


c 


Zr 0.963249-j0.0361306 (yaw tail incidence) L 1-19) 
The time response plots of the yaw normal 
acceleration, angular rate and tail Més aañes are presented in 
Figures 3.19 through 3.21. A close observation of the above pole- 
zero and time response plots of the discrete yaw estimator design 


indicates that they are very close to those of the discrete 


classical design found in the previous chapter. 
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Figure 3.18 Pole-Zero Plot: Uncoupled Yaw Channel Auto rg 
Estimator Design; Discrete Closed Loop Systeme 
Elliptical Airframe 
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с. Simplified Design 


The discrete estimator designed autopilot of the 


Em ueHus section can be simplified by reducing the returning gain 


SO ps. 


appropriate elements of the estimator gain vector: 


resulting simplified estimator yaw autopilot, 


E Ehrough 3.25, 


This can be accomplished by placing zeros into 


-0.0034 
0 
0 
=0:0 17 
0 
-0.0017 
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2. Uncoupled Roll Channel for Cirene 
a. Estimator Gain Vector 
Following the same procedure as in the yaw channel 
case apart from the use of: 
(1) The transpose of the discrete roll system plant matrix 
(2) S-Plane poles slightly faster than those of the continuous 


open loop roll system, Ena 
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b. Design Approach and Analysis 
The discrete estimator designed roll autopilot a n oy 


found by introducing the estimator gain vector of (TIT C SF 


into the origina mo octem: 


114 





1.4 


0.9 


0.8 


0.7 


YAW NORMAL ACCELERATION,HY 
0.5 


0.0 


Figure 3.19 


0.6 





0 60 0 00 00 0 20 00 è e p e è 2 è 29 $$ è e è 2 2 è 2 e e 0 è 2 9 0 è è SS + <- 0 20 è + Sf è Se 2 e e 0 © 1 0 1 OCA A qi. no ooo W WA W nn NW AN W W W c 


® 0000000 ce 000 0g 0° è 00 e è 0 0 e de è a de e e e a KONM NO KON WO GON KON CCC CO COCO O ACO que 0 0 0-0 00000000 = Ө э” ч 


© 00 0 00 00 0 00 00 È e 0 000 000 00 00 00 10 0 1 eh eet eee eee 00 0 0 00 000 0 0 000 ooo oa... o. = = ө ө » 
М < » 


© 00000 000 00 00 0 sh e e e e£ e 6 66 6 1000 0 100010 0000 0 0 0 00 0 00 0000 01 0 0000 1 0 È 0 0 0 20 00 20 211100000 





“een 0 0 0 06 è è è È dee ese de 0 2 è è 0 0 0 00 00 0 0 0 e 0 0 0 0 è 0 È 0 0 00 00 04 0 00 0 0 00 0 0 0000 0 0 0 è 0 È - 0 00000 no 0 0 00 no nO 0 00000000 


1.0 2.0 3.0 4.0 
TIME - SEC 


Yaw Normal Acceleration vs Time, Uncoupled Yaw 


ATA Propor: Erercijmakor Design; Discrete 
EM Eon system, tala ptacal Airframe 


Jecke 


YAW ANGULAR RATE,R 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 


0.05 


0.00 


-0.05 


sep ee aeeeesreeees deg ereggeeeegeee 


See egeep eeeeeeeeepeeee ep ees eepeeeggeeeeesegeeseseeseseeeeee sep ee eeeeeeeeeeeeeegpéë e œ 666 «66 66 66 66666 668 666666 666 66 666 v 


ee eege ege eeeeeeeeeeegggeeëepeeegeeeeeeeeeeeeeeeeeeeggeëepegggeeeeeeeeeeepeeeegeeeeegeegeegeg ees së eeeseeeeseeeeeepeeepee 


e, EE co Co do dc 0 CGC 0 non 0 0 6/5 400 good eloa o 222/2200 e 2/902 c2 ene 





0.0 


Figure 3.20 


1.0 2.0 3.0 4.0 
TIME - SEC 


Yaw Angular Rate vs Time, Uncoupled Yaw 
Channel Autopilot; Beta#mator Реб еп; 
Discrete Closed Loop Verem ELp ЫШ 
Airframe 


(Ab 


0.00 








E : : 
ZA oreorenonnconorncrnanacnarncnnnos piece Lecca enne) Pei sense acess 
! : : : 
D : : : 
Q = å A 
o 5 È : 
SA pie anana een kaaa 
d : 
e : : : 
e : So esis vices skis ; . 
Ali RR ina. AA baa 
: : 
EM MM e aaaea AA i 
? m | | : 
LJ [1] = o 
Hd] Y 
Si: ` on 3 
pa PW... pem EE EEE EEE EE 
i at e : 
D 
o m : : 
CE ee WE lire e |... e AAA 
Gi! DO : 
D : 
s © ab 19 : 
= A EE A ш .:....................... jo eae шш а 
O { T) ! 
4 w T m 
EMEN EE << eco MEE BENED erre + Vesey even e 
= m T + : 
= LL 
A o... 1 c e tn A 
e г 
> | : 
LU 
- 
Seelen eee ee sees ESTEE ETES Een parer erne ai E 
l LI 
a 
- E. KA e Eus 
o 
I 
LU (T) 
3 
© ү 9° Sree eee eee eee Oe ee ere ee eee ee eee eee ee ee eee eT eee ee eee Tee eee eee eee Te ee eee ee ee ee ee ee ee ee ee oq ee eee ee ee ee ee ee ca... 
i s 
94 
` pi LEGEND 
= O_YAW TAIL INCIDENCE dY 
7 
z 
© 
{ 
vo 
a 
0.0 1.0 2.0 3.0 4.0 
TIME - SEC 
Figure 3.21 Yaw Tail Incidence vs Time; Uncoupled Yaw 


Channel Autopilot; Estimator Design; 
Disenete Closed LoepeSystem; Elliptical 
Airframe 


ag 


-1.00 -0.75 


Pi gøre 


-0.50 -0.25 


0 00 0.25 0.50 0.75 


Pole-Zero Plot; Uncoupled Yaw Channel 
Autopilot; Simplified Alterna tor Design 


Discrete Closed Loop System, Elliptical 
Airframe 


118 





1.00 


1.4 


0.5 0.6 0.7 0.8 0.9 


YAW NORMAL ACCELERATION,HY 


0.3 





Fue sac 3 





. 
* 
. 
së eegeeeeeeee WS seeeeegeegegegeegeegeggeeegggeeeeeeeegegeegeg eg eeeeeeeeeeeeeteeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeteeeeeeeeeeeegeegegeeeeeeeeeegeeeeeeeee- 


. 
° 
MAER 
. 


e 
2900 Teeeensgesgeeseeegsssseeeesë se eeëeeeeeeeeeeeeeseëeeeeeeeeeeeeeeseeeeesgëeeegeegeeeeéëgeesepeeeeeegeeeeeeeeee ee eh ee gege eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 
$ 


1.0 2.0 3.0 4.0 
TIME - SEC 


Yaw Normal Acceleration vs Time; Uncoupled Yaw 
Channel 940 Copaslotye Simplaafied Estimator Design; 
Pagka NE sten: Elliptical Airframe 


119 


YAW ANGULAR RATE,R 


0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 


0.15 


oag o 080 80 NO 600 80 00 KG O KO GO 


e Äeeeeeeg se sg eeeeeeeeeeees seess 


Figure 3.24 


oon 00 ho no 80 op n 9 DOG NG BG HON BG 


00000 fee oo 00 05040m BG GG 


e 00° jd 08 0 20 00 2000000 MA4 oo oo 


sos oso ho 00 209000702 00007 dmo 20 





. e 
» 
. e 
. 
e - 
. . 
* e 
* e 
* . 
. . 
. . 
"Ooo B08 000002090000 900089506 0800080600 000 800000 0000 80292 806989820808 040208809050 eg 66666 666 66666 66666 6 6686666 46 6680 eê 66 £ eê 686 666 
. . 
. e 
. . 
* * 
H D 
* . 
. * 
. . 
* * 
. . 
. e 
. . 
eg ee 
. * 
. . 
. * 
* . 
. . 
. . 
e . 
LI . 
* LI 
* . 
* 8 
. . 
16.00000090000000000000000000000000600000000000909000.0...00 + I DI 0 29 2 È 6 66 66 6 66606 +6 60 60 66 666 066 66 0 6669 0 68660 68 66 606 66 +6666 66 64 œ 
. * 
. . 
* . 
. ° 
* * 
. . 
. . 
. 
. . 
* . 
. * 
* . 
Sono oq tosco co O d +00 dd 00 0 LI 0 0 + e «608 6 46 66 66 646 666 < 66 66 66 66 6 e 6 66 6 66 660 66 < Aegeeggeegggeeeeggeeeeggegeeeeeeegeeee 
. * 
. . 
. . 
* . 
. 
. 
. 
e 
* 
. 
. . 
Seege sg gess ges eege 2 e 0 e 0 20 00 029 20 05 299890999090040 00090 0002020900009 0902090090 08999900400 000 8208 280440998 2020099200 2054054250800 
D 
. 
* 
* 
. 
* 
. 
e 
. 
e 
* 
20900000 Ree seg ee geese ee g eg e er eR RR өөө өөө өөө ооо “os. 00 00 00000 CCO 2 0 L8 0 2 8 #6 6 0 è è fg 060 0 a deo 0 è ce dae a dda a + 0 0 cod 0 060 0 0 004 
D 
. 
D 
. 
* 
. 
. 
. 
. 
* 
D 
Lå 
? eee ? 
. 
. 
. 
. 
* 
. 
. 
D 
. 
D 
. 
ec» ee ELLE LL oq coco colsoccooo coros css co soca cosas asd soro 
. 
e 
. 
. 
* 
D 
. 
1 
* 
. 
* 
@eeeeeee 1004001000000 00400 2000000045000 00 cc ceo lo0000000 00000000 c 000000000 1000000011004 
. 
e 
* 
* 
. 
D 
. 
. 
D 
eeeeeee Beesseeseseeeeseeeeeeeeeeeeseseseeeeseleeeeeseeseeeeseegeszegeessssegeeeesssslsssessesgssseseggessegeggggepgesegëgëesg 
. 
D 
ê eee 6 


LEGEND 


e . 
oba? 000000000 057000000090090 00000 9098 20000 000000 000 000500 00 do 
. = 





1.0 2.0 3.0 4.0 
TIME — SEC 


Yaw Angular Rate vs Time; Uncoupled Yaw 
Channel Autopilot: Simple ESU ARR 
Design; Discrete Closed Loop System; 
Elliptical Airframe 


120 


0.00 
4 


«6666666666668 68 68 6 66666 s 668 6 + a 00 dd DÌ da è dI 2 0 00 20 00 0 0008 è 10 00 0 2 do » 1a 020 0 è è 2 20 20 LI 0 66 ê+ êe 666 66686 6 6666 6 e 6666 66 68 66 666 66666 66 666 666 e 6666 66 «686 6 e 00 0005 


-0.01 


è ON GF 6668809 90 9080 0 00 00 00 0 è da è 0000 8084060 069060680. a800 00 000 000 000080000000 900 d0 200000 00 800 04 080200290200000 2000200 00008 29098 809 2988000 


-0.02 


9€9909090909»0900900000090090900090949»9900900V$00209092»»09»09000000909»9900990909»200900200917»220909099090900290»»090900»9009209450090»90200» 80YO PHT G GD 
. 


-0.03 


PAPAA AAA eee sela e ne ao [oo a naje ee dje aje ae ANONS KG (jONAOAM 00 KO SON ais elvis lo ia 6 0 9 o pelo 90 9/9 Vo Do o o 0 ja o jo 0 0/0 0/0 mn 0040 RRR 0/0 0/0 0 00/00 8) 0/0/n/D pio na ere 0/n Dio Bjo o ooo ERER REE 


-0.08 -0.07 -0.06 -0.05 -0.04 


YAW TAIL INCIDENCE,6Y 


-0.09 


ss leese deeeeegeesgeees Pfe ees eegseeeezéësgeeg eg eeeeees eseesgeeesgeseéegessgsgesgëgëegeeegsgeeeegseeegeeseeesëseeessseeesgsseegssesessssssgssgseesgsesssséëëese e 


-0.10 


© 0 ufo v «Kde «66 en «8 ees «e ehe 666 ieee eg 66 66666 66666 66 6 66 66 6 6 +R 6 6 66 66 6 66 66 66 6866 6 6666 66 A 666 66666666 668 66 66666 666 66 eR 666 66 66 6 66 66 6 66 66 6 66 66 6 66 66 6 66 66666 s 


-0.11 


-0.12 





O 
ng 
> 
= 
= 
= 
= 
E 


CIDENCE,6Y | 





-0.13 


e9499 9090600090049 »9092909900090»90900000vo t1, | 0° 00 00 00 000 0000000 


-0.14 


-0.15 


0.0 1.0 2.0 3.0 4.0 


TIME - SEC 


Figure 3.25 Yaw Tail Incidence vs Time; Uncoupled Yaw 
Cosme! suite pilots Simplified Estimator 
Design; Discrete Closed Loop System; 
Elliptical Airframe 


Me 


The pole-zero plot of Figure 3.26 indicates that the 
discrete closed loop system is stable, Since the.z-plane closed 


loop poles were found to be: 
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The time response plots of the roll angle, ange mm 
rate and tail incidence are presented in Figures 3.27 through 
3.29. A close observation of the above pole-zero and time 
response plots for the discrete roll estimator design indicates 
that they are very close to those of the discrete classical desen 
found in the previous chapter. 
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IV. MODERN CONTROL DESIGN MD SN SNE 
COUPLED PITCH AND ROLL CHAN LOR 
CIRCULAR AIRFRAME 
A. GENERAL 
The present chapter deals with the modern control design and 

robustness analysis of the discrete coupled pitch and roll 
channel autopilot for the circular airframe configuration. The 
continuous open loop coupled autopilot whose plant system and 
input matrices are presented in Appendix G is obtained by 
coupling the linear uncoupled pitch [Ref. 10] and roll (Table IV) 
channels. Then, utilizing analog-to-digital conversion by the 
aid of ORACLS program and for a sample period of 0.0125 sem 
the seventeenth-order discrete coupled system with matrices shown 
in Appendix H is formulated. Next, introducing the control-law 
and estimator designs were obtained and analyzed in terms of their 
transient responses and the application of the POPLAR design 
program [Ref. 7]. The POPLAR program is applied in order to 
employ singular value analysis and the use of an optimization 
routine to aid in pole placement control design of the abe 
discussed linear multivariable systems. The robustness of the 
system is also considered by establishing singular value levels 
which correspond to multiloop gain and phase margins determined 
from the universal gain phase system diagram developed by Newsom 


and Mukhapadhyay at NASA Langley. 
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Design Approach and Analysis 


The diserete coupled state-feedback designed autopilot is 


formulated by introducing into the original control system of 
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The pole-zero plot of Figure 4.1 indicates that the 


closed loop poles are: 
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Figure 4.1 Pole-Zero Plot; Coupled Pitch and Roll Channel 
Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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=, = 0.0993086 GRECI 42) 


EM = 0.836555+50.306107 (IV.B.1-13) 
Zz = 0.836553-j0.306107 (IB. 151 4) 
Zig = 0.967788+j0.0334452 mus. 1-15) 
Z< = 0.967788-j0.0334452 (IV.B.1-16) 
Zig = 0.89573+j0. 938133 (CHE. 117) 
Z,3 = 0.89573-j0.938133 (ТОБ. 18) 


The time response plots of the roll angle, angular rate 
and tail incidence are presented in Figures 4.2 through 4.4. A 
close observation of these plots indicate that they are identical 
with those of the discrete uncoupled state-feedback design of the 
previous chapter. 

E»  Robusiness Analysis 

Executing the POPLAR design program of [Ref. 7] with 
inputs the data presented in Appendix I, the minimum additive 
input (MIN ADD IN SV) and output (SVADMO) singular values were 
computed from a frequency range from 0 to 200 rad/sec. 

Figures 4.5 and 4.6 which are plots of SVADMO and MIN ADD 
IN SV versus — indicate that the discrete coupled state- 
pecdoaek design is robust. It is noted that for very low 
frequencies the values of SVADMO are above 0.8680. 

Pay, in terms omwen MIZ3UtyYon results, the ordered 


computed eigenvalues are: 


= 0.07717 (IV.B.2-1) 
= Ey: 09855 (IV.B.2-2) 
2E 0 25923 (IV.B.2-3) 
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Figure 4.3 Roll Angular Rate vs Time; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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Figure 4.4 Roll Tail Incidence vs Time; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circula Airframe 
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Figure 4.5 SVADMO vs Frequency; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; 
БЕ се шше Loop system; Circular Airframe 
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Figure 4.6 MIN ADD IN SV vs Frequency; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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=, = 0.069508 (IV.B.2-4) 
EE 0.8370f-30.3014 (TUNE 2-5 ) 
Se - 0.83707430. 3014 CNB. 2-6) 
Z, = 0.89844 (IV.B.2-7) 
Z, = 0.90064 (IV.B.2-8) 
=, = 0.43996 (IV.B.2-9) 
EN 0. 46207 (IV.B.2-10) 
Z, = 0.96562-j0.03401 (IV.B.2-11) 
Zu = 0.96502+j0.03401 (IV.B.2-12) 
EE 0.99559-50.08307 (IV.B.2-13) 
Z,4 = 0.99559+j0.08307 (IV.B.2-14) 
sii 0.99820 (ERE 2-15) 
a= 1 (IV.B.2-16) 
EL 1.06586 (IV.B.2-17) 


EM GRE TE COUPLED ESTIMATOR DESIGN 
imesh an Approach and Analysis 
Following the same procedure as in the case of the state- 
feedback but for the coupled estimator gain vector the pole-zero 
and time response plots of the coupled estimator design are 
obtained. The time responses are again identical with those of 
the discrete uncoupled estimator design. 
2. Robustness Analysis 
Executing the POPLAR design program of (Ref. 71], Figures 
4.7 and 4.8 are obtained which prove the robustness of the 


system. It is noted that for very low frequencies the values of 


141 







SVADMO are above 0.83976. Finally, in terms of esang in 
results the ordered computed eigenvalues are the same as in the 


coupled state-feedback design. 


142 


0.99976 1.01976 


97976 


SVADMO 
0.91976 0.93976 0.95976 0. 


0.89976 


0.87976 


LEGEND 
SVADMO 





0.839760.05976 


0 10 2 #30 #40 50 6 70 80 90 100 
FREQUENCY 


Figure 4.7 SVADMO vs Frequency; Coupled Pitch and Roll 
Gammel Auntopilot;usBstimator Design; Discrete 
Siroged kecp oysten; Circular Airframe 
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Figure 4.8 MIN ADD IN SV vs Frequency, Coupled Pitch and Roll 
Channel Autoepilot, Esramaveor ean, DUCES 
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A. 


VONC C O AND RECOMMENDATIONS 


CONCLUSIONS 


The goal of the present thesis was the design and analysis of 


weer ete lateral autopilots for application to BTT missiles. The 


following are the principal conclusions based on this work. 


E 


BA 


The continuous and discrete classical designed autopilots 
were proved to have identical performances for the two 
lateral channels. 


The state-feedback and estimator autopilots were introduced 
as additional dynamic designs in order to implement control 
to the original system. Both designs, analyzed in terms of 
their transient responses, were found to meet the desired 
requirements. 


The simplified state-feedback and estimator designs reduced 
some of the returning gain loops, making the system 
simpler, without any significant effects on the system's 
performance. 


The performance of the coupled pitch and roll channel 
autopilot was found to be satisfactory and the overall 
system proved to be robust. 


RECOMMENDATIONS 


iMiworder to improve the simplicity of the overall system, 


more returning gain loops of the state-feedback and estimator 


design must be eliminated. A further investigation then must be 


conducted in order to examine if the performance of the resulting 


design remains unchanged. 
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AP PE Do 


DESIGN REQUIREMENTS FOR UNCOUPEED AUTOrP Rae. 


The requirements for the classical design method of the 
uncoupled channel autopilots [Ref. 6] are the following: 
1. High Frequency Attentuation in Actuator Command Brame 
a. Uncoupled Yaw Channel 
It must be) 15 db at 100 rad/sec and zero angles.” 
attack and sideslip. This requirement will provide sufficient 
high frequency attenuation for 2 30 Hz actuator and for body 
bending modes when high frequency filters are added, but it 
limits the ability of the yaw autopilot to minimize side c ENE 
angle. 
b. Uncoupled Roll Channel 
It must be 2 15 db at 100 rad/sec and zero angle-of- 
attack. This will provide sufficient high frequency attenuation 
for 3 30 Hz actuator and for elastic modes when high frequeme 
filters are added, but tnis requirement limits the speed of roll 
angle response. 
2. Relative Stability for BotneskLateralachemels 
Gain margins >6 db, phase margins }, 30° with a Soa 
12 db апа Ба 
3. Acceleration Time Response 
a. Uncoupled Yaw Channel 


(1) 63% time constant of approximately 0.4 seconds. 
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Em overshnoot ««105. 
(3) Steady-state error need not be zero. 
Uncoupled Roll Channel 

(1) 63% time constant of 0.5 seconds. 
(2) Overshoot « 10%. 


(3) Zero steady-state roll angle error. 
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APPENDIX B 


AERODYNAMIC DATA 


The overall classical design developed by Arrow [Ref. 6] was 
performed for the selected flight condition of Mach number 3.95 
at 60000 feet altitude. The corresponding aerodynamic data 
presented below were taken or derived from the ICAO standard 


atmosphere tables. 


TABLE VII 


RERODINAMIC DATA (M=3.95, BERE 


Temperature, T (OR) 389.988 
Sonie Velocity, a (ft/sec) 968.47 
Pressure, p (1b/ft%) 149.78 
Density, (lb=sec*/ft) 0.0002238 
Velocity, V (ft/sec) 3825.4565 
Dynamic Pressure, q (1b/ft^) 1637.145 
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APPENDIX C 


LINEARIZED ABRODINAMIC DERIVATIVES 


The linearized aerodynamic derivatives at the selected flignt 
condition and about a zero trim angle-of-attack are provided 


below for both airframe configurations [Ref. 6]. 


ТАРЕ Vee a 


LINEARIZED AERODYNAMIC DERIVATIVES 
(M=3.95, H=60kft; a=0°) 


CIRCULAR Bibbie AE 
EA -0.065 -0.043 
EP | -0.025 0.024 
В 
On den 0.016 
Ge. 
EE 20-059 -0.042 
SE 
la 0.031 0.023 
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APPENDIX D 


MISSILE SIZING ANDS MASS be 


In order to provide a realistic missile based on the 
aerodynmamically tested configuration concepts [Ref. 5], the 
models were assumed to be 1/6-scale and the mass properties were 
developed corresponding to mass distribution which might be 
expected for missiles of this size. All the geometric and mass 


properties are presented in the following table. 


TABLE IX 


GEOMETRIC AND MASS PROPERTIES OF MISSILE CONFIGURATIONS 


Greeular Ellipticas 
Length, 1 (in) 168 168 
Max. Diameter (in) 24 
Max. Major Axis (in) 41.57 
Max. Minor Axis (in) 13.86 
e.g. distance BA. (in) TOO OCÓN) 100.8(0.6 ) 
Reference Length, d (ft) 2 2 
Reference Area, S (ft) TI т 
Weight, W (1b) 2525 2475 
Lie (slug-ft^) 4 0 40 
apo res) 810 853 
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АБЕКЕ UD E. 


PROGRAM LOGIC FOR APPLICATION OF ACKERMANN'S FORMULA 


The program logic for computing the control-law gain vector F 
or the transpose of the estimator gain vector KZ, taken from 


КООШО], is given in the following table. 


TABLE 


EROGRAM LOGIC FOR APPLICATION OF ACKERMANN'S FORMULA 


|. Read in t. T.T. and N.: the number of states. 
2. Comment: first we will read in the desired pole locations in the v-plane. convert them 
to z-plane poly nomial coefficients. and construct aid. 

a MES identitv matrix, .N, x A. 

3. ALPHA —] 

5. A — ] 

DE NS roto sep I8. 

EK. Read in pole location 4 us a — yb. 

8. If h = 0. go to step 14. 

p d —lexplalicos bT 

10. As Apa) 

ll. ALPHA — ALPHA ^ (bp: d - ad - A.D 

E A=? 

13. Go to step 6. 

14. A, — expaul) 
IS. ALPHA — ALPHA > (P-A? D 
16. A — A - | 
17. Oo to step 6. 

I8. Comment: now we construct the controllability matrix. 
19. C —1 
20. E—T 
2]. A — ] 
ER- 1I SN BO to step 28. 
E Comment: replace column 4 of C by E: 
24. C[ : A] — E 
25. ATA + | 
26. E — wp - E 
E Go to step 22. 
28. Comment: now solve for the control law. first form ¿7 us the last row of I. 
29. E — HN: | 
AQ. Solve BC = E for B. 
3l. K - B » ALPHA 
EN) 
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M AND INPUT MATRICES OF CONTINUOUS COUPLED PITCH 


AND ROLL CHANNEL AUTOPILOT 
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